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INTRODUCTION 


The  off-line  programming  of  a  two-axis  filament  winding  machine  consists,  in  principle,  of 
the  following  four  general  steps: 

1.  Establishment  of  the  mandrel  geometry  by  the  profile  or  radius  function  r(x), 
where  x  measures  distance  along  the  axis  of  the  surface  of  revolution. 

2.  Establishment  of  the  nominal  winding  path  or  quasi-geodesic  by  the  polar  radius 
function  r0(x). 

3.  Computation  of  T  as  a  single-valued,  periodic  function  of  R,  where  T  is  the  x 
position  of  the  filament  delivery  point  and  R  is  the  aggregate  rotational 
displacement  of  the  filament  delivery  point  around  the  axis  of  the  mandrel. 
Actually,  the  mandrel  rotates  and  the  filament  delivery  point  merely  translates 
back  and  forth  parallel  to  the  mandrel  axis,  but  it  seems  more  useful, 
conceptually,  to  think  of  the  mandrel  as  being  fixed  in  space  and  the  filament 
delivery  point  as  translating  ( T)  and  rotating  (R). 

4.  Proper  scaling  of  the  (R,  T)  winding  data  so  as  to  produce  a  perfectly  uniform 
wrap  which  properly  covers  the  surface. 

This  report  is  mainly  concerned  with  the  proper  accomplishment  of  step  2  so  that  no 
problems  will  arise  with  the  accomplishment  of  step  3.  What  problems  may  arise  while  doing 
step  3?  Well,  if  r0  is  not  defined  properly,  the  (R  T)  winding  data  may  exhibit  loops,  cusps,  or  at 
least  nonunique  T  values  for  a  given  R  value.  That  is,  T  will  not  be  a  single-valued  function  of 
R.  Even  convex  surfaces  are  not  immune  to  this  effect. 


MONOTONE  WINDING  CONDITION 

Step  3  is  accomplished  in  practice  by  computing  R  and  T  as  parametric  functions  of  some 
parameter  such  as  angular  displacement  (of  a  point  on  the  quasi-geodesic  path)  8.  We  obtain 
the  winding  function  T=f(R)  indirectly  by  computing  R=g(9)  and  T=h(8).  The  function/ will  be 
single-valued  if  and  only  if  g  is  a  monotone  increasing  function  of  8.  Hence,  we  will  have 
"monotone  winding"  if  dR/dd> 0.  In  what  follows,  a  prime  denotes  differentiation  with  respect  to 
x  and  a  dot  denotes  differentiation  with  respect  to  8. 
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WINDABILITY 


We  refer  to  a  path  on  the  surface  or  the  corresponding  polar  radius  function  ra  as  being 
admissible  if  r=ra  at  the  two  turning  points  and  rQ<r  at  all  points  in  between. 

If 

p=rr"(rl-r2W0{Urn) 


and 


q-rVrXfl+r'2) 

a  band  of  filament  laid  along  an  admissible  path  will  not  experience  bridging  or  lifting  from  the 
surface  when  the  filament  is  under  tension  provided p> 0  for  all  x.  This  inequality  is  derived  in 
References  1  and  2.  Note  that  an  admissible  path  on  a  convex  surface  can  never  experience 
bridging. 

If  the  inequality  p+q  >0  holds  for  all  x,  the  winding  data  will  embody  a  single-valued 
function  (f ).  Now,  either,  both,  or  neither  of  these  inequalities  may  hold  at  any  given  point.  For 
an  admissible  path  to  be  acceptable,  however,  they  must  both  hold  everywhere.  We  refer  to/»0 
and p+q> 0  as  the  "clinging"  and  "monotone"  inequalities  or  conditions,  respectively. 

Assuming  the  clinging  condition  holds,  the  only  way  in  which  an  admissible  path  may  fail 
to  be  windable  is  when  r'  and  rj  are  of  opposite  sign  and  their  product  is  sufficiently  negative  at 
some  point.  Note  that  at  turning  points  (r=ra),p+q> 0  reduces  to  7+r'ro'>0.  So,  if  clinging 
holds  and  the  product  of  these  derivatives  is  nonnegative,  we  are  theoretically  guaranteed  single¬ 
valued  winding  data.  Also,  if  clinging  does  not  hold  in  some  region  and  the  derivative  product  is 
nonpositive  anywhere  in  that  region,  then  we  are  guaranteed  multi-valued  winding  data. 
However,  if  the  derivative  product  is  positive  throughout  the  region,  we  might  escape  multi¬ 
valued  winding  data  even  though  we  would  still  have  bridging. 


PATH  GENERATION 

In  this  section  we  derive  a  simplified  algorithm  for  generating  quasi-geodesic  paths  on 
surfaces  of  revolution.  The  development  is  more  straightforward  than  that  of  Reference  1. 
Results  of  this  section  are  also  used  in  deriving  the  monotone  winding  condition  in  the  next 
section. 


A  point  on  the  path  on  the  surface  is  given  by  P=ix+jy+kz  where  y=rSinQ,  z=rCos6,  x 
measures  distance  along  the  axis  of  the  surface,  and  6  is  the  angular  displacement  of  P  around 
the  axis. 
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From  y=rSinQ  and  z=rCos6,  we  have 


dy  =drSind+rCosddd 

and 

dz  =drCos8-rSinddd 

If  s  measures  arc  length  along  the  path,  we  have 

ds 2  =dx2+dy2+dz 2  =dx2+dr2+r2dd2 

Now,  if  a)  is  the  angle  that  the  path  makes  with  a  meridian  (0=const),  we  have 

dsSint o  =rdd 

But  the  generalized  Clairaut  relation  defining  a  quasi-geodesic  path  is 

rSinu  =re 

Hence 

dd-—ds 
r 2 

and  therefore 

T 2 

ds 2  =dx2+r'2dx2+—ds2 
_2 


From  which  we  conclude  that 


(  2 

(  2  2\ 

1 

r-r0 

ds2=— 

2 

ra 

r  ~r0 

r2 

A+r,2t 

Mr'2, 

W 


Hence,  we  have  that 

2  =F(x) 


where  dx/d6>0  for  increasing^  and  dx/d6<0  for  decreasing x.  Using  a  three-term  Taylor 
expansion,  we  have 


x=- 


±r 


r-r: 


1+r 


/2 


3 


xto-xoj+Mxe-ej+iMxe-of 


Letting 


A=\+r 


n 


and  differentiating  the  following 

tF(x)  -A  \r*r~2-r2y 


We  have 

*F'(x)=-r'r"A  "2(r*r^-r2y 
+A  *(r 4r  ~2-r 2)  *(2r3r  V  ~2  -r4rj3r'-rr  0 

Since 

x=F(x)F/(x) 

We  finally  have 

x  =-!-[rl(2r2-rZ)-r3rc'1r'-rr/r//(r2-rZ)A  _1] 

& 


Therefore,  we  have  a  simple  second  order  method  for  moving  from  point  to  point  along  a  quasi¬ 
geodesic  path: 
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A=\+r'2 

B  =r  '(2r2-rf)  -r  V e~V '0 
C-rrlr"(r2-r^) 


Ax=xAe+±xA62+0(Ae2) 

2 

Note  that  these  formulas  are  valid  in  the  vicinity  of  turning  points.  This  algorithm  in  conjunction 
with  References  1  and  3  can  be  used  to  generate  winding  data  with  uniform  error. 


MONOTONE  INEQUALITY  DERIVATION 

In  this  section,  we  derive  the  test  condition  which  tells  us  whether  or  not  an  admissible 
polar  radius  function  or  path  is  capable  of  ultimately  producing  single-valued  winding  data. 

Let  P  be  a  point  on  the  path  on  the  surface  and  let  the  taut  filament  be  tangent  to  P  and 
to  the  path.  Let  Q  be  the  other  end  of  the  taut  filament,  called  the  filament  delivery  point, 
residing  some  constant  distance  H  from  the  axis  of  the  surface.  This  configuration  implies  a  two- 
axis  filament  winding  machine,  but  our  resulting  inequality  applies  to  a  three-axis  machine  with 
variable  H  as  well. 

Now, 

Q=P+at 

where  t  is  the  unit  tangent  vector  to  the  path  and  a  is  some  function  of  6.  Since  Q  resides  at 
distance  H  from  the  axis,  we  have 

H2=Q± •  Q± 


where 


v±- component  of  vector  v 
perpendicular  to  axis 


but 
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Therefore 

H2=PX  Px+2aPx-  t+c?tx  tx 
or 


1  = 


r,r,2aP,'h 

H2  H  H 


Knowing  P,  t,  and  H,  we  could  easily  solve  this  equation  for  a  to  obtain  Q,  but  our  intent  here  is 
not  to  actually  compute  winding  data,  but  rather  to  obtain  a  condition  under  which  single-valued 
winding  data  is  guaranteed  to  exist.  So,  if  the  path  can  be  wound  properly  for  a  given  H,  then  it 
must  be  windable  for  any  H\  We  can  therefore  achieve  considerable  mathematical  simplification 
by  simply  letting  H  become  infinite.  From  the  previous  equation,  we  have 

—  =— —  as  H  -  « 

H  «rj| 


where 


llvll  =v/v*v 


Now,  if  X  is  the  angular  displacement  by  which  Q  leads  P,  and  R  is  the  angular  displacement  of 
Q,  we  have 

R=d+X 


and 


PxQx=\Px\\\\Qx\\CosA=rHCosX 


but 

Q^p,+at, 
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Hence, 


P±iPx+at)=rHCos, l 


or 


„  ,  r  a  Px‘  h 

Cos  X= — + - 

H  H  r 


So, 


CosX=Pxtx=4>  os  H  -  oo 


where 


Now, 


dP=P 
ds  s 


but 


P=ix+jy+kz 


Therefore 

Px=jy+kz 

t  -(/y+fc) 


so  that 


p  _  /y+fe  _/y+fe 

Vy2+z2  r 
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We  therefore  have  that 


<f>=CosA=P±- 


-  _  yy+zz 
t>jy2+z2 


But  since 

y=rSin9,  z-rCos6 
y=rSind+rCos6 
z=fCosd-rSin6 


it  easily  follows  that 


yy+zz=rf 


y2+z2=r2+r 


2 


and  that 


Now,  since  we  wish  to  find  a  condition 


= 

\jr2+r2 

equivalent  to  dR/dO> 0,  we  will  need  dk/dO.  Since 
4>=CosX 


we  have 


Differentiating 


<p  =  -Sin  X  X 


— *- 


-<j>\lr2+r2  _-r<j> 
r  r<p 


<f>=r(f2+r2)  2 
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and  since  the  denominator  of  this  fraction  is  always  positive,  we  need 

iz(2  r,2-rrlh)+r2-rr'x  >0 

Now,  recalling  from  the  last  section  that 
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where 


A=l+r 


/2 


B  -r'(2r2-r^)-r3r~lrl 
C=rr'r"(r2-r%) 


we  have 


'  2_  2N 


1^2-  (2r/2-rr/Vr2-^LTB-7)  >  0 

A  rlA\  A> 


Multiplying  this  inequality  by 


A2'! 


gives  us 


but 


A(r2-r£)(2 r/2-rr")+A2r2-r/(AB-C)  >  0 
AB-C=A[r'(2r2-rZ)-r3r;1ryrr,r"(r2-rZ) 


Therefore 

A(r2-rf)(2r,2-rr/f)+A2rl 
-r  'A(j  x(2  r 2  -r^)  -r 3  r  **  rj) 
+rr'2r"(r2-r%)  >  0 

Grouping  terms  involving  r"  and  factoring  out  A  yields 

-rr  "(j  2  -r %)(A  -r'2) 

+A[2rr 2(r2-r02)+Ar2-r/2(2r2-r02)+r3r/r;V^]  >  0 
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which  simplifies  to 


-rr"(r2-r%) 

+A[r/2(2r2-2r2+r2-2r2+r^)+r2+r3r/r~Vj]  >  0 


immediately  yielding 

rr//(r02-r2)+(l+r/2)(rfl2+r3r/r;Vj)  >  0 

as  the  monotone  winding  inequality.  Global  satisfaction  of  this  condition  guarantees  the 
existence  of  single-valued  winding  data. 


WINDING  DATA  GENERATION 


In  this  section  we  derive  the  equations  for  generating  the  nominal  winding  data  that  a 
two-axis  winding  machine  would  need  to  actually  wrap  a  mandrel.  The  development  here  is  not 
new,  but  we  include  it  for  the  sake  of  completeness,  following  the  notation  developed  thus  far. 


Recalling  that 

a2tx  t±+2aP±-  tx+r2-H2=0 


we  solve  for  a,  getting 


h'h  N 


<P:h 


\2 


vV'x 


Hr-r2 

V'x 


>  o 


but 


P=ix+jy+kz, 


s 


Therefore 


p  .t  ~yy+zi-rf 
1  x  s  s 


_y2+z2  _f2+r2 

*-L  .2  .7 

s  s 
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but 


and  letting 


we  have 


So, 


rH 

r 

2 


Pxh_KrH 


K'h 


Kr3r 

a- - + 


N 


K2r6f2  Kr\H2-r2) 
rl 


Kr 


-rr+A 


• o  \ 


2.2  H2-r2 
rzr*+ - 


Letting 


P  =rr,  y= 


H2-r2 
K  ’ 


and 


we  have 


u=-p+\/p2+Y= - - — r 

P+if+Y 
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Recalling  the  general  definition  of  A,  we  have 


H  rH 


H 


Also, 


(XT  X 

T=x+-=x+—?-=x+Kux 


We  can  therefore  summarize  the  equations  for  generating  the  nominal  winding  data. 


•  •  / 

r=xr' 


*=- 


f2+r2 


P=rf,  r  = 


H2-r2 


u  =  ~P+)/p2+Y  if  P  <  0 


u  = - ^ -  if  p  >  0 

P+\/P2+y 


X  -Cos 


r+Kuf 


V  H 
R=0+X 
T  =x+Kux 


SCALING  OF  WINDING  DATA 

In  order  to  produce  uniform  spacing  between  windings  and  ensure  complete  covering  of 
the  surface,  it  is  necessary  to  lengthen  the  nominal  range  of  rotation  for  one  circuit  slightly  and 
to  compute  the  number  of  circuits  and  revolutions  needed  to  cover  the  surface.  We  will  try  to 
do  this  in  a  more  straightforward  manner  than  that  of  Reference  1. 
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If  the  width  of  the  band  of  filaments  is  b,  the  circumferential  cover  afforded  by  the  band 
is  given  approximately  by 


b=cCosa>=c 


N 


Now,  a  circuit  is  defined  as  the  action  of  winding  from  one  polar  parallel  to  the  other  polar 
parallel  and  back  again.  So,  if  we  wind  n  circuits  to  get  two  full  layers  of  coverage,  we  want 

— *27ir  for  all  x 
2 


This  is  equivalent  to 

4lT  L 2  2 
m—Jr  -r0 

But  we  want  this  inequality  to  hold  for  all  x,  so  we  must  have 

4  nJM 

n* — 


where 


M  -Max  r2-r02 


Hence,  we  define  the  nominal  minimum  number  of  circuits  necessary  for  complete  double  layer 
coverage  as 


n.  =ceilind 


(  4  njM 


This  number  of  circuits  used  with  our  nominal  winding  data  has  no  chance  of  producing  a 
uniform  wrap,  however. 

Suppose  the  net  angular  displacement  of  the  filament  delivery  point,  or  the  point  on  the 
path,  or  the  mandrel  itself  for  a  complete  circuit  is  Ra.  Let  i,j,  and  k  be  integers  here.  We  have 

Rg=2ni+p 


where 


Now  define 


^■K) 


where  R,  is  slightly  greater  than  Ra 

R\*Ro 

implies 

■  kp 

J  2  7t 

Therefore  define 

'-“MU) 

It  is  necessary  that  j  and  k  be  relatively  prime,  i.e.,  GCF(j,k)~  1.  If  this  is  not  the  case  for  a 
given  k,  we  simply  do  not  use  that  value  of  k.  The  k  value  is  the  number  of  winding  groups  that 
appear  on  the  surface  and  finally  coalesce  as  a  layer  is  completed.  If  we  scaled  the  nominal 
winding  data  using  Rn  our  winding  pattern  would  close  after  only  k  circuits,  so  we  define 

=/?!+* 


We  want  e  to  be  defined  such  that  after  we  wind  n  circuits,  we  will  have  returned  to  the  very 
same  point  on  the  surface  at  which  we  started,  completing  two  layers  of  uniform  wrap.  That  is 
we  want 


nR2=nRl+ne= 2n 


+ne=2nm 


where  m  is  the  total  number  of  revolutions  of  the  surface  needed  to  wind  a  double  layer.  Hence, 
we  want 


to  be  an  integer.  Moreover,  we  want  GCF(m,n)  —  1.  If  we  wind  k  circuits,  our  net  angular 
displacement  will  be 


kR^  =2  jr(ki+f)+ke 


and  the  filament  will  begin  to  lay  down  alongside  previous  windings.  Hence, 

6  =ke 
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will  be  the  angular  displacement  between  adjacent  windings  in  a  single  layer.  But,  irrespective  of 
the  j  and  k  values,  we  must  have 


if  the  winding  pattern  is  to  close  exactly,  for  the  first  time,  after  exactly  n  circuits  have  been 
wound,  completing  a  double  layer.  Therefore, 

ne _nke _  nd  _  An  _2 
2n  Ink  Ink  Ink  k 


and  we  have 


m  -ni+ 


nj+2 

k 


Now,  if  m  is  to  be  an  integer,  k  must  divide  nj+2  exactly.  Is  there  any  guarantee  that  we  can 
find  an  n  >  =nD  such  that  this  quotient  is  an  integer?  A  fundamental  theorem  of  number  theory 
states  that  given  positive  integers  j  and  k,  there  exist  other  integers  u  and  v  such  that 
uj +vk= GCF (j,k).  In  our  case,  GCF(j,k)= 1.  Therefore,  we  are  guaranteed  integers  u  and  v  such 
that 

^tl=2v 

k 


So,  if  we  add  1  to  -2m,  we  have 

(-2u+l)j+2_w=2v+jl 
k  k 

and  if  we  pick  1  to  be  a  sufficiently  large  multiple  of  k,  w  will  be  a  positive  integer.  Therefore, 
there  are  an  infinite  number  of  values  of  n  >  =n0  such  that  k  divides  nj+2.  We  only  need  the 
smallest  one.  The  reason  that  we  have  defined  n  as  the  number  of  circuits  necessary  to  lay  down 
a  double  layer  is  that  we  want  the  windings  of  the  second  layer  to  be  staggered  exactly  half  a 
bandwidth  from  the  windings  of  the  first  layer  in  order  to  most  efficiently  cover  the  joins  between 
adjacent  windings  of  the  first  layer.  This  staggering  of  the  windings  of  alternate  layers  by 
winding  through  exactly  2vm/n  radians  per  circuit  should  maximize  ultimate  strength.  Our  goal 
then  is  to  find  values  of  j,  k,  m,  and  n  such  that  the  surface  is  completely,  uniformly,  and 
efficiently  covered  for  minimal  m. 
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Now,  the  angular  displacement  between  adjacent  winding  groups  will  be  about  Itt/Ic 
initially,  so  we  will  need 


2jt  ^  .  n 

—  >6  or  k  <  — 
k  2 


Therefore,  in  the  algorithm  to  follow,  we  will  require  that 


s  floor 


The  larger  the  bandwidth,  the  smaller  nc  will  be.  Hence,  no  acceptable  k  may  exist  without 
lowering  the  bandwidth.  This  situation  is  rather  unlikely  to  occur,  however.  The  following 
algorithm  is  suggested. 

”*-skmBn0(i+ 1),  n:=  m 
for  2  a  k 

ifj  >  0  A  /  <  k  A  GCFQJc)  - 1 
then  find  smallest  N  z  ng  such  that  k  divides  Nj+2 

k 


if  GCF(MJf)  =  1  A  [M  <  m  V  (M  =  m  A  N  <  n)] 
then  n:=N,  m:=M  etc. 

When  this  algorithm  terminates,  we  can  define  the  scaling  factor  S  by 

2  nm 


S  =- 


nR„ 


where 


K  ~  Rfinai  R initial 


The  scaled  R  values  are  therefore  given  by 

*s  =  Rinitial+S(R-Riritial) 
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MODIFIED  POLAR  RADIUS  FUNCTION 

The  larger  the  bandwidth  b  is  relative  to  the  mandrel  diameter,  the  more  important  it 
becomes  to  raise  the  polar  radius  function  slightly  thereby  bringing  the  turning  points  in  slight  y 
from  their  nominal  end  positions.  This  is  necessary  because,  ideally,  we  want  to  generate  a  path 
for  the  midline  of  the  band. 

Let  5  denote  meridional  length  and  a  dot  denote  differentiation  with  respect  to  5  here. 
For  a  meridian,  we  have 

ds2  =dx2+dr2  =dx2+r'(x)2dx2 


Hence 

x  =(\+r'(xf)  2 
jc  =  -r/(x)r//(jc)(l+r/(x)2)-2 

A  three-term  Taylor  expansion  for  x  as  a  function  of  s  is  given  by 

x(s)  -  x(fl)+x(a)(s-fl)+jx(a)(s-fl)2 

=  x(a)  +(1  +r/(x(a))2)'2(5-a)-ir'(x(a))r//(x(a))(1  V^a))2)' V*)2 


Let  the  length  of  the  mandrel  be  L  and  let  M  be  the  total  length  of  the  meridian  between 
nominal  turning  points  atx=0  andx=L.  The  new  left  turning  point  (s=b/2,  a=0,x(a)-U)  will 

then  be 


*, = j=(i  +rWv\-\rW\m  V(0)2r2-j 


b  b2r/(0)r"(0) 
2v/l+r/(0)2  8(l+r/(°)2)2 


and  the  new  right  turning  point  ( s=M-b/2 ,  a=M,  x(a)-L)  will  be 

*r.  £.<1  >r%L?Y^ 

b  b2r'(L)r"(L) 

2/r^w  8(i+rW 
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These  approximations  to  the  new  turning  points  should  be  adequate  for  most  cases.  Now  we 
must  raise  the  polar  radius  function  so  that  it  agrees  with  the  radius  function  at  the  new  turning 
points. 

Define 

6,=r(x}-r0(x^  Sr=r(x)-r0(x) 


<$(*)=  <5, 


(x-x)(6r-6l) 


xr-x, 


If  we  denote  the  modified  polar  radius  function  by  rx,  we  then  have 

rjx)=  rg(x)+6(x) 


I  I  &T~&1 

rUx>  ro  (*)  + - 


xrxi 


and 

We  would  therefore  use  r ^  in  place  of  ra  in  the  equations  of  all  sections  previous  to  this  one. 
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